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Abstract 

We discuss the problem of static chaos in spin glasses. In the case 
of magnetic field perturbations, we propose a scaling theory for the 
spin-glass phase. Using the mean-field approach we argue that some 
pure states are suppressed by the magnetic field and their free energy 
cost is determined by the finite-temperature fixed point exponents. In 
this framework, numerical results suggest that mean-field chaos expo- 
nents are probably exact in finite dimensions. If we use the droplet 
approach, numerical results suggest that the zero-temperature fixed 
point exponent 9 is very close to &dL m In both approaches d = 3 is 
the lower critical dimension in agreement with recent numerical sim- 
ulations. 
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1 Introduction 



One of the most interesting open problems in spin glasses regards a correct 
understanding of the nature of the low temperature phase, i.e. the spin-glass 
phase n, 0]. Spin glasses are characterized by a strong freezing at a certain 
critical temperature. Below that temperature a complete description of the 
nature of the static phase is still missing. 

During the last years there have been several developments in the field, 
the most well known being the mean- field theory 0. Unfortunately mean- 
field theory has revealed a complex theoretical structure which is very ob- 
scure when applied to non exactly solvable models for which some kind of 
perturbation theory is needed. 

There are other approaches to spin-glasses which are known as phe- 
nomenological droplet models, a complete description of them has been given 
by D. S Fisher and D. Huse |7|, |(|. The main idea underlying these models is 
that the spin-glass behaviour is governed by the zero-temperature fixed point 
in the renormalization group equations. [J. Up to now it seems that the 
Parisi solution to mean-field theory is essentially correct. It has passed the 
stability analysis || and gives also a correct description of the thermodynam- 
ics, in agreement with the numerical simulations. It is not clear what is the 
correct description of the spin-glass phase in short-ranged models. Droplet 
models are expected to be a good description of the low temperature phase 
mainly in the case of low dimensions. But droplet models are not suited to 
describe the physics of high dimensional systems and particularly mean-field 
theory. 

The complexity of the replica approach is found when studying the spec- 
trum of fluctuations around the Parisi solution. The full set of gaussian 
propagators has revealed a very complex structure H and the obtention of 
the one-loop corrections to the mean-field equations makes progress slow. 
The main difficulty of this task is the enormous number of sectors within 
replica space which contribute to the one-loop correction. This explains also 
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why finite-size corrections to the main thermodynamic functions are still un- 
known in mean- field theory JTIJ. To all these problems should be added also 
the fact that, up to now, the major part of the computations have been done 
only close to T c within the Parisi approximation. 

In this work we will try to introduce a different approach to the problem 
which can help in understanding the nature of the spin-glass phase. The 
main idea of the approach is to try to look for one order parameter whose 
spectrum of fluctuations is easier to take into account. The chaos problem 
was proposed some time ago and concerns the chaotic nature of the spin 



glass phase flT2"| , |TT|| . The term chaotic can be misleading since it can evoke 
different meanings. In this context, we prefer to use the word static chaos. 
By this we mean that a small perturbation of the Hamiltonian is enough to 
reshuffle the Boltzmann weights of the different equilibrium configurations. 
One constructs a system which is the sum of two Hamiltonians, the initial plus 
the perturbed one. The full system lives in a larger phase space and allows 
for a new order parameter. This order parameter is the overlap between the 
equilibrium configurations of the initial system and the perturbed one. This 
new order parameter has a longitudinal spectrum of fluctuations without 
zero modes and hence is stable. The associated correlation functions to this 
order parameter decay exponentially to zero with a caracteristic correlation 
length. 

The nature of the chaos problem is also interesting concerning numerical 
techniques like simulated annealing where the the change of the temperature 
has to be considered as a perturbation to the system. In this case one wants 
to reach the ground state after a progressive cooling of the system. Let us 
suppose that the spin glass behaves chaotically against temperature changes. 
Then, the equilibrium configurations should reorganize completely for any 
small change of the temperature and a slow cooling would be useless. A 
small change of temperature would be considered like a new quenching and 
the system would be always strongly far from equilibrium. Fortunately, as 
we will discuss later for this particular problem, if there is chaos against 
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temperature changes then it is small and finite-size corrections ensure that a 
high degree of correlation between the equilibrium configurations at the two 
temperatures is preserved. 

The work is divided as follows. In the following section we will present 
a quantitative definition for chaoticity and we will introduce different type 
of chaotic perturbations. We will also present predictions from mean-field 
theory and phenomenological droplet models. Section 3 is devoted to the 
study of a particular perturbation, i.e. chaoticity against changes of magnetic 
field. Starting from the mean-field approach we propose a scaling behaviour 
within the spin glass phase. We discuss also our predictions in the framework 
of the droplet approach. Section 4 is devoted to numerical simulation results 
and finally we present our conclusions and a discussion of the results. 

2 A working definition for chaos 

The idea underlying the chaoticity of the spin-glass phase relies on the fact 
that it is a marginal phase |9[|. The fact that the spin-glass phase is not fully 
stable means that a small addition of energy to the system is able to change 
completely the statistical weights of the equilibrium states with a very small 
cost of free energy (of order 1/N compared to the supplied energy to the 
system, where N is the size of the system). 

Marginality is one of the outstanding results in mean-field theory of spin 
glasses. It is also a feature of phenomenological droplet models and in general 
it is related to the fact that in the spin-glass phase spatial (time) correlation 
functions decay very slowly with distance (time). This decay is not far from 
a power law in the most general case. The full reorganization of equilibrium 
states in spin glasses after a small perturbation is a natural feature in mean- 
field theory. In this case there is an infinite number equilibrium states and 
all of them contribute to the partition function but with a different weight 
[fnjl . This is because they have equal free energies per site except differences 
of order 1/N. Any small but finite addition of energy to the system is enough 
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to redistribute these small free energy differences reshuffling the weights of 
the different pure states. 

In droplet models there exists the concept of overlap length (it is some- 
times denoted Lat or L^h according to the case if the perturbation is a 
change of the temperature or the magnetic field). Droplet models suppose 
that there is only one equilibrium state. When the system is perturbed, the 
correlation functions reorganize completely in a scale of distances larger than 
the characteristic overlap length. It is clear anyway that the overlap length 
in these models has to be always smaller than the correlation length and 
coincide only in the limit in which the perturbation vanishes. 

In what follows < ... > and (.) mean thermal and disorder average respec- 
tively. Now we want to give an appropiate definition of what is static chaos. 
For simplicity we will consider Ising spin glasses even though the definition 
can be generalized to other models. Let us suppose an Ising spin glass system 
with Hamiltonian Hi [a] . Then we apply a perturbation P to the system and 
the new Hamiltonian for a different copy of spins {rj} is given by 



We consider now a full Hamiltonian which is the initial system Hi[a] 
plus the perturbed one H 2 [t], i.e. H[a,r] = Hi[a] + H 2 [t]. The phase 
space has been enlarged and we can consider a new order parameter which, 
for example, in the case of Ising spin glasses, is given by the overlap (<7jTj) 
between the equilibrium configurations of the system Hi with the equilibrium 
configurations of the perturbed one H 2 . When there is no perturbation, i.e. 
P = 0, this is the usual order parameter of the spin glass with Hamiltonian 
H = Hi + if 2 . 

We now define the chaoticity parameter r by: 



where (cf of ) denotes the order parameter evaluated taking two copies a, b 
of the unperturbed system Hi and similarly for (rf rf) of the perturbed sys- 



H 2 [r}=Hi[r]+P[r] 



(1) 



r(P) 




(2) 
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tern B.2- The thermal average of the order parameter in the numerator is 
performed with the full Hamiltonian H. In principle, this order parameter is 
equal to one if the perturbation is zero. This is trivial because r is the order 
parameter of the spin glass normalized to itself. Chaoticity in spin glasses 
reflects the fact that any small but finite perturbation P causes the parame- 
ter r to fall abruptly to a value smaller than r = 1. Obviously this can only 
happen in the thermodynamic limit because, for a finite size N, the chaos 
parameter r will always be a smooth function of the perturbation P. This 
means that one has to perform the thermodynamic limit before applying the 
perturbation P. More precisely, the spin-glass phase is chaotic if 

lim lim r(P) < 1 (3) 

It is also possible to define the adimensional quantity 

a = j3L (4) 

where the numerator is obtained by averaging over the Hamiltonian H[ar] 
and the denominator is the order parameter for two copies a and b of the 
same system The difference between the adimensional quantities a 

and r is only an appropiate normalization. In fact, the definition of chaos 
given above in eq.(Q) also holds in case of the parameter a. The necessity 
to distinguish among the parameter a and the parameter r is important 
for certain types of perturbations. For example, in the case of temperature 
changes the order parameter (erf a!-) is very sensitive to the temperature and 
vanishes at the critical point. Let us suppose the initial Hamiltonian is in 
the low temperature phase and we change the temperature by putting the 
system close to the critical point. The chaos parameter a vanishes because the 
numerator in eq.([|) vanishes close to T c and the denominator remains finite. 
On the contrary, the chaos parameter r of eq.(|2[) has in the denominator a 
term which also vanishes at T c and normalizes appropiately the numerator. 
Because r measures correctly the overlap among equilibrium configurations 
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it is the appropiate parameter to deal with in case of temperature changes. 
The difference between the chaos parameters a and r is not important in the 
case of magnetic field changes and other types of perturbations. 

Let us now discuss what happens in the case of ordered systems. As an 
example we take the standard Ising model in a finite number of dimensions. 
Let us suppose that we are in the low temperature phase, at temperature 
T below the critical point, and let us take as a perturbation a small change 
of the temperature. At the temperature T the system has a spontaneous 
magnetization m. There is only one equilibrium configuration with a fraction 
m of the spins pointing in a certain direction. When we change a little bit 
the temperature by a small quantity AT, the mean number of spins which 
point in that direction (i.e. the magnetization) changes linearly with AT at 
least for AT small. In this case one sees inmediately that eq.(^) gives the 
value r(AT) = 1 for any small change of the temperature and the system is 
not chaotic. We represent an equilibrium state by an N dimensional vector 
v = (mi, m2, •••) Wat) where the are the local spin magnetizations. In the 
Ising model a slight change of the temperature modifies only the length of 
this vector but not its direction. The chaos mechanism in spin glasses is 
driven by the fact that as soon as we perturb the system this vector v suffers 
a sudden rotation because the weight of the different states are changed. In 
most cases, any small perturbation makes the vector v to become orthogonal 
to its previous value and r(P) vanishes for a finite perturbation P. 

There are many examples of perturbations that one can apply to the sys- 
tem. As was mentioned in the introduction, one can change the temperature 
or change the magnetic field. These are among the most studied pertur- 
bations in the literature. But one can imagine other kinds of perturbation 
like for instance changing the realization of disorder. In this finite 
fraction of the couplings is changed (for instance, in case of symmetric 
distribution of couplings, this change could consist in reversing the sign of 
the perturbed couplings J^). One can also imagine to add a small ferromag- 
netic or antiferromagnetic part to the couplings. In these cases, spin glasses 
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seem to behave chaotically against these perturbations. As an example, we 
show in figure 1 how the parameter r decreases with the size N, for the 
case of the SK model for two different perturbations. The first perturbation 
corresponds to a very large change of temperature (AT = 0.4) for an initial 
temperature T = 0.5 = T c /2 (the perturbed system is at T = 0.9 which is 
very close to T c but always in the spin-glass phase.) The other perturbation 
is the application of a small magnetic field h = 0.2 to a system initially at 
zero magnetic field and T = 0.6 (the AT line lies at h ~ 0.4). The results 
for N less than 20 have been obtained by calculating exactly the partition 
function, the remaining ones using Monte Carlo simulations. From figure 1, 
the system seems much more sensitive against magnetic field perturbations 
than temperature changes. This is clear also if we observe that under the 
temperature perturbation, even though it is very strong because it puts the 
perturbed system close to the paramagnetic phase, the equilibrium configu- 
rations at both temperatures still retain a high degree of coherence (r ~ 0.7). 

In the following, we will focus on the study of a particular perturbation, 
which has turned out manageable in order to understand its effects in the 
spin glass phase: the case in which the perturbation consists in applying a 
small magnetic field to a spin-glass at zero field. This has been the subject 
of previous research, specially by I. Kondor in the case of mean-field theory 
[ p~5|j . One could also study the case in which the system is at a finite field 
in the spin-glass phase and the field is slightly changed. This problem is 
more subtle than the previous one in which the system is initially at zero 
magnetic field. The main reason is that (at least for short ranged systems) 
we do not know if the spin-glass phase survives to a magnetic field. If the 
spin-glass phase survives to the magnetic field then we expect (as predicted in 
the mean-field approach) that chaoticity will be present in a magnetic field. 
In the other case (and this is the prediction of droplet models), the system 
would then be always in the paramagnetic phase and chaos should not be 
present. Then, according to eq.(Q), r(P) = 1 + O(Ah) would be continuous 
for Ah = 0. 
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In the case the system is initially at zero magnetic field mean-field ap- 
proach and droplet models agree in that they both predict that the spin- 
glass phase is chaotic. More specifically r(h) (we use the intensity of the 
applied magnetic field h for the pertubation P) is zero for any finite h. 
But the main mechanism which makes the spin-glass phase chaotic is very 
different in both pictures. In phenomenological droplet models the spin- 
glass phase is marginal: the correlation functions decay very slowly with 
the distance and the correlation length associated to the two point function 
C(x) = (cr(O)cr(x)) 2 is infinite. When a magnetic field is applied the spin- 
glass phase is destroyed and the correlation length becomes finite. It is given 
by 0: 

t~(q E Ah 2 )^ (5) 

with qEA the Edwards- Anderson order parameter and 9 the thermal exponent 
which gives the characteristic energy scale If of droplet excitations of typical 
size L. This means that all excitations of droplets of sizes larger than a 
certain length £ will be suppressed by the field. The exponent 9 is a zero 
temperature exponent (it is determined by the zero temperature fixed point of 
the renormalization group equations) and it is expected to be constant in the 
low temperature phase. In the critical point the associated thermal exponent 
9 C is determined by the finite temperature fixed point of the renormalization 
group equations and is related to the critical exponents by 9 C = 2b|±zz where 
1] is the anomalous dimension exponent an d is the dimension (even though 
it has been argued that at low dimensions there appears a new exponent 9 C 
[|T^1). In general, we expect that 9 C is smaller than 9 above the lower critical 
dimension and both vanish in the lower critical dimension. 

Mean-field theory approach gives a completely different mechanism of 
chaoticity. After applying a small magnetic field, the spin-glass phase is not 
destroyed. We suppose that the effect of the magnetic field is the dissapare- 
ance of a large number (infinite) of equilibrium states. This mechanism is 
easy to visualize by taking into account the correct order parameter for spin 
glasses which is the distribution P(q). Its physical meaning was explained 
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some time ago [[18], [19] and it gives the probability density that two pure 
states a and (3 have a common overlap q a p = q. This common overlap corre- 
sponds to the scalar product of the local spin magnetization in both states. 
At zero magnetic field the function P(q) is symmetrically distributed around 
q = and non zero within the interval (— q max , q m ax)- in a magnetic field the 
reversal symmetry a — > —o is broken and P(q) is non zero only for q positive 
and larger than a minimum value q m in- The value of q ma x is nearly indepen- 
dent of the magnetic field (this approximation, which works extremelly well 
close to T c , is called the Parisi- Toulouse hypothesis |T6[). In some sense the 
effect of the magnetic field is to suppress those equilibrium states a which 
had overlaps q a p with the other remaining states (3 smaller than q m i n . Within 
the usual picture of the spin-glass phase in mean-field theory JTJJ there is an 
infinity of states with a few number of them dominating the Gibbs measure. 
This infinity of states lay in the tips of an ultrametric tree and the effect of 
the magnetic field corresponds to progressively cutting those branches which 
generate the states which are suppressed. The suppression of the states also 
conserves the ultrametricity property. The understanding of how pure states 
a are suppressed by the field according to their statistical weight w a is still 
an interesting open problem. 



3 Chaos in magnetic field 

This section is devoted to the study of chaos in case of an Ising spin glass 
initially at zero field after turning on a magnetic field. We are interested in 
the case of a <i-dimensional Ising spin glass with random Jjj couplings defined 
by the Hamiltonian 

H = - JijOiGj - h^Oi (6) 

(m) * 

where the couplings Jy are quenched variables distributed according to a 
probability function P(J) of zero mean and finite variance. The interaction 
is restricted to nearest neighbours and h is the magnetic field. The Ising spins 
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(jj take the two possible values ±1 and live in a <i-dimensional hypercubic 
lattice. In the limit d — > oo one expects to converge to mean-field theory, 
i.e. the SK model. In the SK model all spins interact one to each other and 
the couplings Jjj are normalized by a factor 1 / ^/N where N is the number 
of spins. 

3.1 The case of mean- field theory 

This question was adressed some time ago by I. Kondor |L5]]. Let us consider 
two copies of the same realization of disorder, one at zero magnetic field and 
the other one at finite magnetic field h. The full Hamiltonian of the problem 
is given by: 

H[a,r] = H l [a]+H 1 [r]-hJ2r i (7) 

i 

with 

Hi = - J H °i a i ( 8 ) 

l<i<j<N 

This problem can be directly solved using the standard replica trick for the 
full Hamiltonian J2a=i, n H a where H a is given by eq.(^) and a is the replica 
index which runs from 1 to the full number of replicas n (at last one takes 
the limit n — > 0). Now one applies the replica trick 

b£^=]im^— ^, (9) 

which yields the expression 

T} = J dP ab dQ ab dR ab exp(-NA[PQR}) (10) 

with 

A[PQR] = ^(P* b + Ql + 2Rl b ) + J2Rl-logTr aT e W ((3 2 J2Pa b <Ta<T b 

a<b a, a <b 

+ P 2 Y,QabTaT b + (? £ R ab O a T b + 1 £ R 

a<b a^b a 

+0/»E 7 J ( n ) 
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In this way, one is able to reduce the problem in terms of a lagrangian 
A[PQ R] with three order parameters corresponding to the different overlaps 
among the two copies, i.e P ab = (a a a b ); P aa = , Q ab = (r a r 6 ); Q aa = and 
Rab — (&a T b)- F° r finite h there is an inmediate solution for the equations 
of motion which is given by the P ab and Q ab Parisi matrices with zero and 
magnetic field h respectively and R ab = 0. The free energy of the whole 
system is the sum of the free energy of one copy at zero magnetic field plus 
the free energy of the other copy with field h. This is a solution because it 
gives the full free energy of two uncoupled systems. The order parameter 
associated to R is 



In order to study the stability of this solution one computes the spectrum of 
fluctuations. The full set of fluctuations is very complex. For instance, within 
the subspaces generated by the diagonal subblocks P and Q, it corresponds 
to the spin-glass spectrum derived by C. De Dominicis and I. Kondor [|]. 
Only the fluctuations around R = (the off-diagonal subblock) are those 
which are physically relevant to the problem because they measure spatial 
correlations between states corresponding to the two Hamiltonians, the initial 
and the perturbed one. In mean-field theory there are no distances and we 
want to obtain the spatial behaviour of the system within the mean-field 
approximation. This can be done using a Ginzburg-Landau approximation 
by introducing spatially dependent order parameters in the effective action 
of eq. fliTf) . Now the order parameters P,Q,R depend on the space variable 
x and we add to the action A a kinetic term of the type YL a <b 9 fgz • The 
spectrum of fluctuations is contained in the the momentum space propagator. 
This is given by the Fourier transform G(p) of the correlation function 



where (•) means averaging over disorder and ((•)) is the usual thermal average 
over the Hamiltonian eq.(0). 



(12) 




(13) 
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The problem of computing the propagator reduces to the diagonalization 
of a hierarchical matrix of the Parisi type. The full expression has been 



reported in |Tj|. Its singular part is given by 

rq max rQma* p 2 + 1 + A(g)A(Q) 

G(p)= i dq L, dQ (f + l-X( q )X(QW (14) 



with 



PQmax 

X(q) = (3(1 - q max + / dqx(q)) (15) 

J a 



where f3 is the inverse of the temperature. The same expression applies in 
the case of X(Q). Here q(x) and Q(x) are the order parameter functions for 
the spin glass at zero and h field respectively. 

The correlation function eq.(^) decays to zero for large distances x with 
a characteristic length £ which is given by the minimum eigenvalue of the 
stability matrix. This eigenvalue is non-zero for finite h which demonstrates 
the stability of the R — solution. The correlation length £ (which is the 
inverse square root of the minimum eigenvalue) diverges like (1 — \(Q m in)) 

2 2 

Close to T c = 1 we have Q m in ~ h~s. This gives £ ~ h~3 which diverges when 
h -> 0. 

The stability of the R = solution implies that the system is chaotic. 
This means that r(P) (as given by eq.(|2])) always vanishes for finite h like 
jj, where iV is the size of the system. The result that £ diverges when h 
goes to zero is rather natural because in that case the perturbation vanishes 
and the two copies are identical. Then the correlation length £ is the spin- 
glass correlation length which is infinite because there is marginal stability. 
In some sense, there is a first order phase transition at h = where the 
probability distribution associated to the order parameter q defined in eq. ( |T2] ) 
changes from a delta function peaked in q = to the usual order parameter 
distribution for the spin glass pD| . Also the critical point is chaotic but in 
this case the correlation lenth diverges like £ ~ h~^. We remind the reader 
that in the paramagnetic phase there could not be chaos of the type defined 
in eq.(H). This is because the correlation length £ would never diverge but 
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only converge smoothly to its corresponding finite value at zero magnetic 
field. 

Now we turn to the behaviour of the propagator G(p) of eq.(|T4|) in the 
limit p — > 0. Using the known expressions [^T|] for q(x) and Q(x) close to T c 
in eq. (|I3D we obtain a divergent expression for G(0). Its most divergent part 
is 

G(0) ~ / ^ r (16) 

jQmin (1 - A(Q))§ 

which gives G(0) ~ p~ A ~ £ 4 ~ /i~f . This is not new and this result can 
also be obtained from the study of the intravalley gaussian propagators as 



derived in |22[. We can define a certain kind of non-linear susceptibility by: 

= = G(o)~/r* (17) 



This susceptibility can also be written \ni = N (q 2 ) with q given in eq. fljjD . 
Using the fact that R = is a stable solution altogether with eq.flTTD, the 
following scaling behaviour holds 



a = r ~ 



f(Nhz) (18) 



This result will be derived in the following section using scaling arguments 
and will be also generalized to short-range models. 



3.2 Scaling theory of chaos with magnetic field 

Next we want to give a precise physical meaning to the correlation length 
£. As commented in the previous section, the spin glass phase is marginal 
with an infinity of equilibrium states, none of them having a characteristic 
correlation length. Under a small magnetic field, a lot of states are suppressed 
and the correlation length £ is finite. We interpret £ as the new typical 
correlation length of the states which have been suppressed. This is the 
natural continuation of what comes out in the critical point. In this case there 
is only one marginal state. After applying a magnetic field, the correlation 
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length becomes finite and the system goes into the paramagnetic phase. In 
the spin-glass phase the suppressed states acquire finite correlation length 
and their free energy increases respectively to the remaining ones. We are 
still within the spin-glass phase because the remaining states dominate the 
partition function and they still have infinite correlation length. In the spin- 
glass phase all equilibrium states are non equivalent. Some of them have a 
much higher statistical weight. This means that only those states {a} which 
give overlaps {q af 3 < q m i n V/3} are simply erased by the magnetic field. When 
all states are suppressed we reach the AT line ||13|| . This can only happen 
in case when an infinity of equilibrium states coexist in the low temperature 
phase. 

If we want to be more precise we have to generalize these ideas to the case 
of short-range models. Two basic assumptions are enough to this aim. The 
first one concerns the physical interpretation on the effect of the magnetic 
field on the equilibrium states. The second one uses information in the 
critical point to understand what happens in the spin-glass phase. More 
precisely, we suppose that the cost in free energy of the dissapearing states 
in the spin-glass phase, scales in the same way as in the critical point. We 
argue that the low temperature spin-glass phase is determined by the finite 
temperature fixed point of the renormalization group equations. This is the 
contrary assertion of droplet models in which the spin-glass beahviour is 
governed by the zero-temperature fixed point. Our assumptions give exact 
results in mean-field theory. The existence of some critical properties in the 
low temperature phase has been also seen in a different context. For example, 
it has been proved that the exponent which characterizes the decay of the 
tail of the P(q) around q = q ma x freezes below the critical point in mean-field 
theory fnj| . In short-range Ising spin glasses there are also numerical results 
which suggest that the freezing of some critical exponents really takes place 



in the low temperature phase [23. 24 . 

At the critical point we know that the singular part of the free energy 
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(per site) is given by 

fsing ~ h 2 q ~ gM (19) 

Here g is the usual order parameter defined in eq.fll^D, d is the dimension 
and [q] is the dimension of the operator Q a b in units of the inverse of the 
correlation length, the value [q] is connected to the critical exponents /3, v 
and rj via the relation [g] = f = d ~1 +v ■ 

Now we generalize this expression to the case in which replica symmetry 
is broken, i.e. in the spin-glass phase. First of all, we need a general ex- 
pression for the singular part of the free energy which is invariant under the 
permutation group of the different replicas. The most easy expression of this 
type is 

fsing ~ E Qfl (20) 
a<b 

where the exponent [q] is given by the critical exponents. 

One can easily derive the correct behaviour of the singular part of the 
free energy for an order parameter q(x) of the type shown in figure 2. To 
obtain the correct singular part of the free energy corresponding to the states 
which are suppressed by the field we have to take the difference of eq.( 
with h ^ and h = 



fsing = EQa b (h)-Y,Qa b (0) 
a<b a<b 



| d d d d 



+1 



(Qh ( X ) 1o ( X )) d-X — Xmin Qmin ~ Qmin (21) 

The main ingredient that we have used in this derivation is the fact that 
the order parameter q(x) in short-range models is characterized by a con- 
tinuous part plus a plateau. Under the application of a magnetic field a 
new plateau appears with g(x) = q min and x min ~ q m %n- This last result 
is connected with the fact that the order parameter distribution P(q) at 
zero magnetic field is finite for q = 0. Because P{q) = (where x(q) 
is the invers of the q(x)) this means that q(x) ~ x for x close to zero. In 

2 

the critical point the previous derivation applies with x m i n = 1 and q ~ h~s 
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where 5 = fzf^ 2 - For droplet models the same derivation is valid but now 
%min — 1 and -P(O) vanishes in the thermodynamic limit like L~ e with 6 the 
zero-temperature fixed point exponent. The singular part of the free energy 
scales like £ being given by eq.([5j). 

Now we apply eq.(pTD to mean-field theory. Mean- field critical exponents 
together with rj = give [q] = 2. Because q min ~ ftJ we obtain f sing ~ /it 
and the global singular free energy scales like N ftJ . Because the parameters 
a and r are adimensional we reproduce the scaling behaviour of eq. (|T8|) . 

3.3 Estimate of the AT line and the lower critical di- 
mension 

The first result which comes out from the previous subsection is that d = 4 
plays a role as a special critical dimension. This deserves some explanation. 
The upper critical dimension in Ising spin glasses is 6. Above 6 dimension 
the critical exponents coincide with the mean-field ones. These exponents 
are associated with the order parameter Q(x) corresponding to the overlap 
a a (x)a b (x) between two copies a, b with the same Hamiltonian. The correla- 
tion length associated to the two point function (Q(0)Q(x)) diverges at the 
critical point and remains infinite in the low temperature phase. The chaos 
correlation length is associated to the two point function (R(0)R(x)) and 
corresponds to a different order parameter R(x) = <j(x)t(x) which couples 
two systems with different Hamiltonians. We argue that the exponent of the 
chaos correlation length £ associated to R(x) lies in a different universality 
class of that to which the order parameter Q(x) belongs. 

We can find the appropiate upper critical dimension associated to the 
criticality of chaos. From eq. fll8D and using £ ~ we obtain an argument 
of the scaling function for a of the form L/£ in four dimensions. This means 
that d u = 4 has the role of an upper critical dimension. 

In the most general case we can introduce the exponent A defined by 
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Qmin ~ h*. We expect the following scaling to be satisfied 

2(d+M) 

a = r ~ f(Nh~^T) (22) 

The value [g] depends on the critical exponents and this scaling contains only 
one non critical parameter (A) and thus is easily measurable in a simulation. 

The exponent A is theoretically unknown and there is no numerical pre- 
diction on its value. From the value of A one obtains the AT line using the 
condition q min ~ q ma x- Because q max ~ t 13 with r = -^p^ the AT line is 
given by the equation h ~ r^~. In mean-field theory (3 = 1 , A = 3 gives 

3 

h ~ T2. In d = 4 depending on the value of A and /3, a different expression 
is found. This and the special case d = 3 are left as a discussion in the next 
section. 

One can also estimate what is the lower critical dimension d%. In fact, we 
expect that the scaling eq. (|22]) for T c = should be of the form a = f(N h 2 ). 
This gives = X[q), i.e. [g](A — 1) = d. The exponent A should diverge as 

d approaches d\ because q is discontinuous at T c = when a magnetic field is 
applied. So [q] = 0, i.e, di — 2 + r] = which is the usual relation determining 
the lower critical dimension p5| (in principle this relation should at least 
be satisfied for Hamiltonians with a countinuous distribution of couplings). 
Furthermore, in case d = di one expects £ ~ h d i . Because £ ~ h 3 for 
d = d u = 4 this means that di = 3 at zero order of approximation or mean- 
field level. We call it mean-field level or zero order because in this case we 
suppose the exponent for the correlation length £ does not vary between 
d u = 4 and d = 3. The fact that 4 and 3 are very close assures that this is a 
good approximation which is probably exact. 

We should now recall that all these predictions have to be apropiately 
modified for droplet models. For these models there is no transition in a 
magnetic field. The condition di — 2 + rj = also applies and the exponent 
6 is well approximated by the result 9 = ^p. This is in agreement with 
the numerical results of the following section. Recent numerical simulations 



for case d = 2 show that the zero temperature exponent 9 ~ —0.46 [26[ is 
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surprisingly close to the chaos prediction —0.5. This suggests that in the 
framework of droplet models also di — 3 and the previous expression for 9 
are probably exact. 

4 Numerical results 

In this section we present Monte Carlo numerical simulations in order to 
test these ideas. We should note that the chaos parameters a and r defined 
in this work and all the scaling laws based on them are computable using 
standard numerical simulations. The standard technique is to consider two 
parallel Monte Carlo simulations, one for the system H\ and the other one 
for the perturbed system H 2 = P[Hi]. The first copy is at zero magnetic 
field while the second one has a magnetic field h. Both copies evolve in time 
and, once they have thermalized, one computes the corresponding order pa- 
rameters. Since one is interested in scalings within the spin-glass phase, the 
main difficulty is that samples have to be equilibrated in the low temperature 
phase where metastability is very strong. All the results in this section are 
for small lattices and we have paid attention that they are fully equilibrated. 
The general schedule of the simulation is as follows. An initial cooling is per- 
formed until the first copy at zero magnetic field thermalizes at the working 
temperature and the second copy thermalizes with an applied magnetic field 
equal to a maximum value h max . Then, the first copy evolves without any 
perturbation and the field of the second copy is progressively decreased step 
by step down to zero. In general, for each different value of the magnetic 
field of the second copy, a long enough thermalization is done after which 
statistics is collected. Then, the order parameters a and r of eqs.([|) and 
eq.(|2|) can be computed. 

The Hamiltonian under study is given by eq.(^]). In all simulations we 
have used the heat-bath algorithm and spins are updated sequentally. In 
case of short-range models we impose periodic boundary conditions. The 
distribution of the coupling J is discrete (the J^- can take the values ±1 with 
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equal probability). If there is a finite temperature phase transition we expect 
universality to apply (anyway see P7| ) and the results for discrete couplings 
should be equivalent to the case in which the distribution is continuous. Our 
main goal is now to test scaling laws of the type eq.(^). Scaling fits work 
well if we use the parameter r or the parameter a (in all cases they differ very 
slightly, approximately by 5 per cent). Then, we will present results only for 
the parameter a. 

We now show the results in case of mean-field theory. The results for 
different magnetic fields ranging from 0.2 up to 1.0 at T = 0.6 are shown in 

3 

figure 3 for several sizes. We show the parameter a versus hN& (we have 
chosen this argument instead of N hi in order to compare directly these 
mean-field results with those corresponding to short-range models.) There 
is an agreement with the prediction of eq. (|T8|) . At the critical point T = 1 
the appropiate scaling argument is N h 3 and in order to compare with the 
scaling law eq.(|T8|) of figure 3 we show results for T c = 1 in figure 4. If in 
figure 3 we plot the chaos parameter a in function of N hi (instead of h ivl) 
we discover that the scaling functions of figures 3 and 4 are clearly different 
suggesting that the criticality of chaos in the critical point and in the low 
temperature phase are in a different universality class. 

Next we present results for the case d = 4. Figure 5 shows the parameter a 

2 

as a function of L hs . This is the mean-field scaling which is in full agreement 
with data. Simulations were performed at T = 1.5 (T c ~ 2.05 pEfl ) which 
is ~ 0.7T C . Metastability effects are very strong and thermalization is more 
difficult (in the sense that one needs more thermalization steps) than in 
the SK model case. Error bars are not shown because they are very small 
(of order of the size of the symbols). The agreement with the theoretical 
prediction is good. 

Then we can derive results for the AT line. Using eq.(^) we get A = 4.2 
which gives q m in ~ h 0A8 . In the critical point taking 77 ~ —0.25 and using the 
hyperscaling relation 5 = ^zfqr where q ~ one gets the result 5 ~ 3.6. In 
the critical point, q scales with the magnetic field with a larger exponent q ~ 
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h ' 55 and one has the impression that this is a general feature at any dimension 
(in the mean-field case q scales linearly with the field in the critical point 

2 

while the minimum overlap scales like hs in the spin-glass phase.) Because 
the critical exponent (3 ~ 0.6 and rj = —0.25 ±0.1 the corresponding AT line 
should scale like h ~ r L 3±o.i w j 1 j c j 1 j s c l se to the mean-field theory result 
(even though there is no reason that it should coincide). Unfortunately we 
have no means to test if this prediction is correct, mainly because the question 
of the existence of the AT line is still unsolved |29], [!(], |3T|. In the framework 
of droplet theory we can derive the correct value of the zero-temperature 
exponent 9. It gives 9 = 0.5 (in the critical point the finite-temperature 
exponent 9 is 9 C ~ 0.43). Figure 5 also gives valuable information in the 
case there is no AT line. The chaos correlation length should correspond to 
the correlation length of the spin-glass in the paramagnetic phase as given 
in eq.(^). In the case of d = 4 we obtain £ ~ 5h~z if we estimate £ as 
the distance over which the chaos paarmeter a decreases by an order of 
magnitude. In order to search numerically for the existence or not of phase 
transition in magnetic field at T = 1.5 one should study lattice sizes L > £ 
where £ is given by the previous expression. 

We analyze now the data for d = 3. Simulations were performed for 
the ±J nearest-neighbour Ising spin glass. Recent numerical simulations 



so large correlation length, we expect that the system will have some kind 
of pseudocritical behavior for small lattices. In fact, standard numerical 
simulations for small sizes show that T c ~ 1.2 with rj ~ — 0.1,f3 ~ 0.5 



thermodynamic limit, simulations for small lattices should be sensitive to 
the pseudocritical behaviour and finite-size scaling for the chaos parameter a 
should give information regarding this pseudocritical point. What comes out 
is very interesting and has been plotted in figure 6. The mean-field result 

2 

£ ~ h~A fits very well the data. This is in agreement with the fact that 
d\ = 3 which is the value for the lower critical dimension if mean-field theory 




But, due to the 



[|33| , [34]] . This means that, even if there is no true phase transition in the 
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is exact. If di = 3 then any finite T belongs to the paramagnetic phase 
which is characterized by the true finite correlation length £y (in order to 
distinguish it from the chaos correlation length £). The chaos correlation 

2 

length £ would behave like /i~s in the regime £ << £t- in the regime £t ~ £ 
the value of the chaos correlation length £ should progressively match the 
value of £y and remaining finite. We can derive in this regime of sizes the 
location of the pseudocritical AT line. From eq.(p2|) we derive q m i n ~ ^°' 26 
which gives h ~ r for the critical exponent (3 ~ 0.5. This result is not far 
from experimental determinations of the AT line in bulk CuMn spin glasses 
|35|j (where a scaling of the form h ~ r 1,8 with an exponent slightly larger 
than the mean- field result 3/2 is compatible with experiments.) 



5 Conclusion 

It seems that the study of static chaos in spin glasses can give interesting 
predictions of the nature of the spin glass phase. The information obtainable 
from the subject is great because of the different ways one can perturb the 
system. In this work we have focused in a magnetic field perturbation. In this 
case it is possible to establish a physical picture in which states are supressed 
by the action of the magnetic field. By 'suppression of the states' we mean 
that these states increase their free energy and do not contribute any more 
to the partition function. Similarly to what happens at the critical point 
where there appears a finite correlation length after applying a magnetic 
field, these suppressed states acquire a finite correlation length (the chaos 
correlation length). To prove this result we should know what is the real 
mechanism of the modification of the free energies of the different equilibrium 
states. This means to understand how states are suppressed by the magnetic 
field depending on their statistical weights. This is an interesting analytical 
problem in mean-field theory which is possibly not out of reach. 

Using this ideas we have been able to derive a scaling behaviour for the 
chaos parameters a and r which depends on the critical exponents at T c . 
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We argue that some critical exponents survive in the spin-glass phase which 
means that the low temperature phase is governed by the finite-temperature 
fixed point. Curiously this is the opposite assertion of droplet models in which 
the spin-glass behaviour is governed by the zero-temperature fixed point. The 
complete understanding of the correct description of the low temperature 
phase in spin-glasses is one of the major still open problems. Mean-field 
approach yields the value 2/3 for the chaos correlation length exponent. We 
expect this exponent to be exact down to d = 4 which is the upper critical 
dimension if hyperscaling applies for the singular part of the free energy of the 
supressed states. Our numerical results are in very good agreement with this 
prediction. Surprisingly this 'mean-field' behaviour seems exact down to d — 
3 which should correspond to the lower critical dimension. In the framework 
of droplet picture our numerical results suggest that the relation 9 = is 
probably exact. Both approaches predict that 3 is the lower critical dimension 
even though we expect some kind of pseudocritical behaviour in the regime 
in which lattice sizes are smaller than the true correlation length. This 
pseudocritical behaviour is expected also with finite magnetic field giving a 
pseudocritical AT line. Our results predict a transition line in agreement with 
some experimental results. For d > di — 3 the physics is determined by the 
existence or not of a spin-glass phase with magnetic field. A definite answer 
on the existence or not of the AT line in finite dimensions is a prioritary task 
in order to clarify the controversy on the real nature of the spin-glass phase 



(see [36] for some recent numerical results.) 

Now we say few words in case of temperature changes. In this case, 
as shown in figure 1, chaos is much weaker. This is a interesting result 
which finds also a natural explanation in the context of mean-field theory. 
This result has already been shown in [|Tj| doing the same kind of analysis 
as has been performed in section 3.1. Namely, when the initial system at 
temperature T and the perturbed one at temperature T' both lie within the 
spin-glass phase then X(q) in eq.flTop is always equal to one which corresponds 
to an infinite correlation length. This means that if chaos exists then it is 
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marginal and the chaos parameter r for any finite perturbation AT goes to 
zero like N~ a with an exponent a < 1. Obviously this result does not exclude 
the possibility that chaos is absent and a = 0. Why chaos is much weaker 
in case of temperature changes than for magnetic field perturbations can be 
intuitively understood if one imagines that by lowering the temperature new 
equilibrium configurations emerge from previous ones and that the system 
suffers a continuous bifurcations into new states. In this degree of 

coherence has to be preserved between the new and the old states and this 
is in agreement with mean-field calculations on the chaos problem. 

It would also be very interesting to understand the chaotic nature of the 
low temperature phase in other spin-glass models, random-field problems 
and the vortex-glass phase in superconductors. The chaos approach could 
reveal as a good starting point to obtain (like happens in Ising spin glasses) 
the lower critical dimension for several models in which there is still much 
controversy. 

From the study of chaoticity in spin glasses we also expect to give some 



hint regarding some real dynamical experiments in spin glasses [37, 38]. Cy- 
cling temperature experiments show that by lowering the temperature some 



degree of correlation is preserved between the probed states|39], pEOfl. Even 



though experimental spin glasses never reach equilibrium we think that a 
correct answer to the statics is relevant to a qualitative understanding of the 



effect of perturbations in the out-off equilibrium relaxations [O, f42|. There 



are also in course some cycling magnetic field experiments which we hope 
will be in agreement with the main conclusions of this work. 



Acknowledgements I gartefully acknowledge very stimulating conversa- 
tions with A. J. Bray, S. Franz, D. Huse, M. A. Moore and G. Parisi. This 
work has been supported by the European Community (Contract B/SC1*/915198). 
I am grateful to J. Kurchan, E. Marinari an M. Picco for a careful reading 
of the manuscript. 



23 



References 

[1] M. Mezard, G. Parisi and M. A. Virasoro in "Spin Glass Theory and 
Beyond" (World Scientific, 1987) 

[2] K. Binder and A. P. Young, Rev. Mod. Phys. 58 (1986) 801 

[3] D. Sherrington and S. Kirkpatrick, Phys. Rev. B 17 (1978) 4384 

[4] W. L. Mc. Millan, J. Phys. C: Solid State Phys. 17 (1984) 3179 

[5] A. J. Bray and M. A. Moore, Phys. Rev. Lett. 58 (1987) 57 

[6] G. J. M. Koper and H. J. Hilhorst, J. Physique (France) 49 (1988) 429 

[7] D. S. Fisher and D. A. Huse, Phys. Rev. B 38 (1988) 386 

[8] C. De Dominicis and I. Kondor, Phys. Rev. B 27 (1983) 606 

[9] C. Dc Dominicis and I. Kondor, J. Physique Lett. 45 (1984) L-205 

[10] G. Parisi, F. Ritort and F. Slanina, J. Phys. A: Math. Gen. 26 (1993) 
247 , 3775 

[11] A. J. Bray and M. A. Moore, J. Phys. C: Solid State Phys. 17 (1987) 
L463 

[12] S. R. McKay, A. N. Berker, S. Kirkpatrick, Phys. Rev. Lett. 48 (1982) 
767 

[13] J. R. de Almeida and D. J. Thouless, J. Phys. A: Math. Gen. 11 (1978) 
983 

[14] M. Mezard, G. Parisi, N. Sourlas, G. Toulouse and M. Virasoro, J. Phys. 
(France) 45 (1984) 843 

[15] I. Kondor, J. Phys. A: Math. Gen. 22 (1989) L163 



24 



[16] G. Parisi and G. Toulouse, J. Phys. Lett. (France) 41 (1980) L361 

[17] M. Nine and H. J. Hilhorst, Phys. Rev. Lett. 68 (1992) 2992 

[18] C. De Dominicis and A. P. Young, J. Phys. A: Math. Gen. 16 (1983) 
2063 

[19] G. Parisi, Phys. Rev. Lett. 50 (1983) 1946 

[20] A. P. Young, Phys. Rev. Lett. 51 (1983) 1206 

[21] G. Parisi, J. Phys. A: Math. Gen. 13 (1980) 1101,1887 

[22] T. Temesvari, I. Kondor, C. De Dominicis, J. Phys. A: Math. Gen. 21 
(1988) L1145 

[23] G. Parisi and F. Ritort, J. Phys. A: Math. Gen. 26 (1993) 6711 

[24] J. C. Ciria,G. Parisi and F. Ritort, J. Phys. A: Math. Gen. 26 (1993) 
6731 

[25] C. De Dominicis, I. Kondor and T. Temesvari, Int. J. Mod. Phys. B7 
(1993) 986 

[26] N. Kawashima, N. Hatano and M. Suzuki, J. Phys. A: Math. Gen. 25 

(1992) 4985 

[27] L. Bernardi and I. A. Campbell, received preprint 

[28] R. R. P. Singh and S. Chakravarty, Phys. Rev. Lett. 57 (1986) 245 

[29] A. J. Bray and S. A. Roberts, J. Phys. C: Solid State Phys. 13 (1980) 
5405 

[30] E. R. Grannan and R. E. Hetzel, Phys. Rev. Lett. 67 (1992) 907 

[31] J. C. Ciria, G. Parisi, F. Ritort, J. J. Ruiz, J. Physique I (France) 3 

(1993) 2207 

25 



[32] E. Marinari, G. Parisi, F. Ritort preprint \cond-mat 9310041 



[33] R. N. Bhatt and A. P. Young, Phys. Rev. B 37 (1988) 5606 
[34] A. T. Ogielsky, Phys. Rev. B 32 (1985) 7384 

[35] G. G. Kenning, D. Chu, B. Alavi, J. M. Hamman and R. Orbach, J. 
Appl. Phys. 69 (1991) 6240 



[36] M. Picco and F. Ritort, prperint \cond-mat 94030711 



[37] Ph. Refregier, E. Vincent, J. M. Hamman and M. Ocio, J. Physique 
(France) 48 (1987) 1533 

[38] L. Sandlund, P. Svendlindh, P. Granberg and P. Norblad, J. Appl. Phys. 
64 (1988) 5616 

[39] J. M. Hamman, M. Ledermann, M. Ocio, R. Orbach and E. Vincent, 
Physica A 185 (1992) 278 

[40] H. Rieger, preprint cond-mat 940241 

[41] J. P. Bouchaud, J. Physique I (France) 2 (1992) 1705 

[42] L. F. Cugliandolo and J. Kurchan, Phys. Rev. Lett 71 (1993) 173 

[43] J. M. Hamman, private communication 



26 



FIGURE CAPTION 



Fig. 1 Chaos parameter r in the SK model for two different pertur- 
bations. In one case the system is at T = 0.6,h = 0. and we apply a field 
h = 0.2, T staying constant. In the other case the system is at T = 0.5 and 
we increase the temperature by AT = 0.4. Error bars in the second case are 
smaller than the size of the symbols. More details in the text. 

Fig. 2 Mean-field theory order parameter function q(x) in the spin- 
glass phase with magnetic field. It is characterized by a minimum overlap 
q m i n , a maximum overlap q max and the corresponding breakpoints x m i n ,x max . 

Fig. 3 Chaos with magnetic field in the SK model at T = 0.6. Field 
values range from /i = 0.2upto/i=1.0 for the smaller sizes and up to h = 0.4 
for the largest ones. The number of samples range from 200 for N = 32 down 
to 40 for iV = 1632. Typical error bars are of order 5 per cent in all cases. 

Fig. 4 Chaos with magnetic field in the SK model at T — T c — 1. Field 
values range from h — 0.1 up to /i ~ 0.8. The number of samples range from 
200 for N = 32 down to 50 for N = 736. 

Fig. 5 Chaos with magnetic field in the 4d ±J Ising spin glass at 
T = 1.5. Magnetic field values range from h — 0.1 up to h — 1. The number 
of samples is approximately 100 for all lattice sizes. Typical error bars in 
this case are of the size of the symbols. 
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Fig. 6 Chaos with magnetic field in the 3d ±J Ising spin glass at 
T = 1.5. Magnetic field values range from h = 0.1 up to h — 1. The number 
of samples is approximately 200 for all lattice sizes except for L = 7 in which 
there are 100 samples. Typical error bars are shown in case L = 6. 
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